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Abstract — Electronics components have and increasingly 
critical role in avionics systems and for the development of 
future aircraft systems. Prognostics of such components is 
becoming a very important research filed as a result of the 
need to provide aircraft systems with system level health 
management. This paper reports on a prognostics 
application for electronics components of avionics systems, 
in particular, its application to the Isolated Gate Bipolar 
Transistor (IGBT). The remaining useful life prediction for 
the IGBT is based on the particle filter framework, 
leveraging data from an accelerated aging tests on IGBTs. 
The accelerated aging test provided thermal-electrical 
overstress by applying thermal cycling to the device. In-situ 
state monitoring, including measurements of the steady- 
state voltages and currents, electrical transients, and thermal 
transients are recorded and used as potential precursors of 
failure. 1,2 
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1. Introduction 

Predicting the remaining useful life (RUL) of electronic 
components like MOSFETs, BJTs, and IGBTs is one of the 
most challenging frontiers of Prognostic Health 
Management (PHM) systems. Such components form the 
backbone of avionics systems that play an ever-increasing 
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critical role in on-board, autonomous functions for vehicle 
controls, communications, navigation, and radar systems. 

During high-voltage operation, large internal electrical 
fields build up within power devices, which can increase the 
rate of degradation. For prognostic purposes the externally 
observable characteristics of the device, such as 
measurements of the collector current/voltage during duty 
cycle on/off states of the IGBT/power-FET that show 
degradation trends towards thermal runaway or latch-up, are 
correlated to physics-of-failure mechanisms like hot carrier 
injection, electro-migration, etc., and used in a particle 
filtering (PF) framework to carry out RUL prediction. 

This prognostic exercise demonstrates the feasibility of 
detecting failure precursors in semiconductor device 
behavior and using them in an intelligent prediction 
framework to derive RUL estimates. The fact that the 
particle filer algorithm provides a RUL probability 
distribution (pdf), instead of a single Mean-Time-Between- 
Failure (MTBF) value, makes the interpretation of the 
prognostic result more useful and intuitive. 


2. Background 

The research in prognostics of IGBTs has been approached 
from different angles. A study of failure precursors for 
prognostics has been presented in [1]. This work focused on 
devices that had been aged under thermal overstress and 
identified precursors of failure that could be measured and 
computed from external signals. An accelerated aging, 
characterization and operational scenario simulation system 
was presented in [2]. This system allows the aging under 
different operational conditions providing thermal and 
electrical overstress to devices, while monitoring in-situ key 
operational parameters that could serve as precursors of 
failure. The development of physics based models of aging 


1 



and fault progression has also been considered. The work 
presented in [3] focuses on modeling the aging effects at a 
system level, by studying electrical power drivers using 
IGBTs. 


Failure Mechanisms Review 

Some major intrinsic faults, relevant to the transistor 
physics, include dielectric breakdown, hot carrier injection, 
and electromigration [4, 5]. Some major extrinsic faults, 
relevant to the transistor packaging, include contact 
migration, wire lift, die solder degradation and package 
delamination [4, 6]. 

Dielectric breakdown occurs when a strong electric field 
induces a current channel through a previously insulated 
medium. Acute dielectric breakdown is typically the result 
of electrostatic discharge (ESD) and junction over-voltage. 
Time Dependent Dielectric Breakdown (TDDB) refers to 
the breakdown of gate oxide caused by chronic defect 
accumulation in the Si02 insulator. TDDB is shown to be 
advanced by increases in electric field strength [5]. A strong 
electric field may impart energy into an electron or a hole 
which then becomes a so “ hot carrier ” due to the high 
kinetic energy [4]. Hot carriers have sufficient energy to 
tunnel and become trapped in gate oxide, and they are the 
primary cause of TDDB and contribute to device failures 
under normal operating conditions. Electromigration is a 
result of high current densities in silicon interconnects 
causing migration of metals. Formation of metal voids on 
interconnects can cause open circuits or high resistive paths, 
which in turn can result in poor performance or circuit 
malfunction [5]. 

Contact migration forms metal voids between external 
contact metals and the silicon. As metal voids grow, the 
aluminum or other metals can diffuse down to the silicon. 
This in turn can cause metal spikes to form deep in the 
silicon region which results in shorting the p-n junctions 
[7]. Wire lift occurs when the bond between the package 
wires connecting to the silicon die fail. Wire lift has been 
identified as a dominant failure mode in high power IGBTs 
[6]. Die solder degradation is another prominent package 
related fault. Solder attaching the silicon die and package 
heatsink develop cracks and voids due to thermal expansion 
mismatch between materials during expansion and 
contraction [6, 8]. 

Accelerated Aging Methodologies 

Thermal stress and electrical stress are the most common 
aging methodologies. Thermal cycling and chronic 
temperature overstress are prevalent thermal stress methods, 
with thermal cycling among the most prevalent accelerated 
aging methodology in electronics. Thermal cycling subjects 
devices to rapid changes in temperature differentials 
causing thermal expansion and contraction. Die solder 
degradation and wire lift are associated strongly with this 


aging method. Thermal overstress, another prevalent 
method, subjects devices to high temperatures for extended 
periods of time. TDDB is accelerated under high 
temperatures [9] and transistors have exhibited temperature 
dependant lifetimes accelerated by this mechanism [10]. 
IGBTs aged with self heating have shown changes in 
current ringing characteristics during switching [3]. 

Electrical overstress can be induced though transient and 
steady-state methods. Transient methods include electro 
static discharge (ESD), inductive switching and 
electromagnetic pulses. ESD is a leading cause of gate 
oxide failure [11] and hard switching of inductive loads, 
causes voltage spikes which can cause significant damage to 
drain-source junctions [12]. Steady-state methods include 
chronic over-voltage and over-current. Applying high gate 
voltages [13], setting gate voltage (V g ) to maximize drain 
current [14], and applying current overstress across the 
drain [15] have been shown to induce hot carrier and TBBD 

[9]. 


3. Remaining Useful Life (RUL) 

Particle filtering is a sequential Monte Carlo (SMC) 
technique for implementing a recursive Bayesian filter using 
Monte Carlo (MC) simulations. It is primarily used for state 
estimation and tracking. The mathematical formulation for 
PF methods have been discussed in [16]. The basic idea is 
to develop a nonparametric representation of the system 
state pdf in the form of a set of particles with associated 
importance weights. The particles are sampled values from 
the unknown state space and the weights are the 
corresponding discrete probability masses. As the filter 
iterates the particles are propagated according to the system 
state transition model, while their weights are updated based 
upon the likelihoods of the measurement given the particle 
values. Resampling of the particle distribution is done when 
needed in order to prevent the degeneracy of the weights. 
For state prediction purposes the same PF framework can be 
used by running only the model-based particle propagation 
step until the predicted state value crosses some 
predetermined end-of-life (EOL) threshold. The predicted 
trajectory of each particle then generates an estimate of 
RUL, which can be combined with the associated weights to 
give the RUL pdf. The process is broken down into an 
offline (learning) and an online (tracking and prediction) 
part. During offline analysis, regression is performed to find 
representative ageing curves. Exponential growth models, 
as shown in Eqn. 1, are then fitted on these curves to 
identify the relevant decay parameters like C and A: 

0=CQxp(-At) (1) 

where 6 is an internal model parameter of interest. More 
details of the PF framework used here can be found in [17]. 

The state and measurement equations that describe the 
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semiconductor aging model are given below: 

Z* = z*.i.exp{-A k (t k -t k . 1 )}+ io k 

A k = A*.i + v k 
[z k 9 AaJ 

y k = z k +v k (2) 

where the vector z, in our case, consists of the off state 
exponential time decay constants for the IGBT collector 
current (I C e), and matrices C and A contain their aging 
decay parameters C and X values respectively. The z and A 
vectors are combined to form the state vector x. The 
measurement vector y comprises the time decay parameters 
inferred from measured data. The time index is denoted by 
k. The values of the C and A vectors learnt from regression 
can used to initialize the particle filter. The noise samples 
v and v are picked from zero mean Gaussian 
distributions whose standard deviations are derived from the 
given training data, thus accommodating for the sources of 
uncertainty in feature extraction, regression modeling and 
measurement. System importance resampling of the 
particles is carried out in each iteration so as to reduce the 
degeneracy of particle weights. This helps in maintaining 
track of the state vector even under the presence of 
disruptive effects like unmodeled operational conditions (in 
our case, high temperatures). 

The system description model developed in the offline 
process is fed into the online process where the particle 
filtering prognosis framework is triggered by a diagnostic 
routine. The algorithm incorporates the model parameter as 
an additional component of the state vector and thus, 
performs parameter identification in parallel with state 
estimation. Predicted values of the time decay parameters 
are compared against end-of-life thresholds to derive 
EOL/RUL estimates. Figure 1 shows a simplified schematic 
of the process described above. 



Figure 1 — Particle Filter Framework 


4. Experiment Details 

The IGBT behaves essentially like a MOSFET and a BJT 
connected in a Darlington configuration. The input of the 
IGBT is modeled by the MOSFET, while the output is 
modeled by a BJT. As a combination of these two devices 
the IGBT takes on their characteristics as well as some 
additional ones caused by their interaction. For example, it 
exhibits a tail current while turning off which is modeled by 


a decaying current flowing through the BJT after the 
MOSFET has been turned off. This is shown in Figure 2. 



Figure 2 - IGBT Voltage and Current Waveforms 


At high temperature thermal cycling these current decay 
curves change due to hot carrier injection. The trends in 
these curve shifts can act as a precursor to catastrophic 
latch-up and subsequent thermal runaway conditions. Figure 
3(a)-(c) shows the shift in the off-state current decay as 
aging progresses, finally ending in a latch-up condition as 
shown in Figure 3(d). 
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Figure 3 - (a), (b) and (c) Shift in IGBT Off-state 
Current Decay with Aging; (d) Latch-up Condition 

Accelerated Aging Settings 

Thermal overstress accelerated aging was conducted on 
IGBTs using the agile accelerated aging, characterization 
and scenario simulation system for gate controlled power 
transistors presented in [2]. An International Rectifier 
IRG4BC30KD IGBT with a 600V/15A current rating in a 
TO220 package was attached to the transistor test board 
with no external heatsink. The collector-emitter junction 
was connected in series with a load power supply and a 
0.2Q load resistor. A 50Q resistor was placed between the 
gate driver and the IGBT gate for current measurement. An 
infrared sensor was attached to the IGBT case for 
temperature measurement. The gate signal was chosen to be 
a Pulse Width Modulated signal with amplitude of 10V, 
frequency of 10 KHz and a duty cycle of 40%. The aging 
process was controlled by a hysteresis temperature 
controller with set points of 268°C and 270°C, switching 
the gate voltage ON and OFF as a control action. The load 
power supply voltage was 10V which resulted in a current 
of 8 A. An additional temperature threshold controller, with 
a set point of 305°C, was programmed to turn off the load 


power supply and end the experiment in the event of 
thermal runaway and latching failures. The IGBT under 
consideration was aged for ~210 minutes until latch-up 
resulting on thermal runaway. 


5. Results 

During the model learning process, we extract the tail 
current sections from the I C e waveform, transform them to 
the log domain and then fit 3 rd order polynomials, as 
expressed below: 

Ice (0 = exp{p/ + P 2 t 2 + P 3 t + P 4 }. (3) 


The equivalent fits in the time domain are denoted by the 
cyan lines in Figure 3(a)-(c). The fit parameters Pi, P 2 , P 3 
and P 4 are then tracked over the entire aging duration. It is 
to be noted that the parameters are highly correlated with 
each other (as denoted by the correlation coefficient matrix 
shown below in Eqn. 4), and consequently, we use only Pi 
as the parameter to track. 


1.0000 -0.9964 0.9873 -0.9747] 

-0.9964 1.0000 -0.9972 0.9899 

0.9873 -0.9972 1.0000 -0.9978 ' 

-0.9747 0.9899 -0.9978 1.0000J 


(4) 


Coeffof Polyfit Term P^t® 



Figure 4 — Trending of Curve Fit Parameter Pi 


Overall, the parameter shows an exponential growth rate 
(negative) as shown in Figure 4, indicating an Arrhenius 
aging process perhaps. Near the end, the curve shows some 
anomalous behavior as the IGBT approaches the latch-up 
condition. Originally the Arrhenius equation described the 
temperature dependence of the rate of a chemical reaction. 
Currently it is perceived largely as a empirical relation 
which can be used to model the temperature variance of 
many thermally-induced processes. 

For running the particle filter framework, we use an 
exponential growth model for Pi like the one shown in 
equation 2. Since we do not have separate learning and 
testing datasets, we do not use the regressed C and X values 
from the model learning process to jump start our PF 
algorithm. Instead, we use 0 and 1 as the initial values for Pi 
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and yl P] . The EOL threshold is chosen as -2.5x10 s , which is 
approximately the value of Pi at the end of its exponential 
behavior. The prediction time is arbitrarily chosen to be 
51.875 minutes. 

Figure 5 shows the performance of the PF in tracking and 
prediction modes. The inset plot shows the 2 Pj parameter 
estimation. After the prediction timeline, the mean of the 
last 100 2 Pj values is used for the propagation of each 
particle according to the state transition model (first line in 
equation 3). No more computed Pi values are used to update 
the particle weights. The propagated particle values are 
compared against the EOL threshold to compute the 
distribution of EOL times. A mixture of Gaussians is then 
fitted in a least square sense to these EOL values to plot the 
magenta EOL pdf. The RUL pdf is simple the prediction 
time subtracted from the EOL pdf. It is to be noted that 
despite the generic initial particle values, modeling 
uncertainties and system noise, the PF performance is very 
good. 

6. Conclusions 

This prognostic exercise demonstrates the feasibility of 
detecting failure precursors in semiconductor device 
behavior and using them in an intelligent prediction 
framework (e.g. particle filtering) to derive RUL estimates. 
The fact that the PF provides a RUL pdf, instead of a single 
Mean-Time-Between-Failure (MTBF) value, makes the 
interpretation of the prognostic result more useful and 
intuitive. 


The results presented here incorporate simple exponential 
growth models to explain the aging behavior of IGBTs. 
Although the prediction results are very good, such an 
approach has limitations in trying to apply the aging 
parameters and EOL thresholds learnt from one device to 
another. In order to achieve that, we need to explicitly 
derive physics-of-failure based aging models, identify their 
parameters from externally observable device characteristics 
and then use those models in a PF framework to carry out 
prognostics. In such a case, aging behavior will be linked to 
the changes in the internal model parameters which can then 
be tied to failure mechanisms like hot carrier injection, 
electro-migration etc. The focus of our future work will be 
to devise aging experiments to fulfill the above objectives. 
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